Planckian Birth of the Quantum de Sitter Universe 



J. Ambj0rn a,c A. G6rlich b J. Jurkiewicz b , and R. Loll c 

a The Niels Bohr Institute, Copenhagen University 
Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark, 
email: ambjorn@nbi.dk 

b Institute of Physics, Jagellonian University, 
Reymonta 4, PL 30-059 Krakow, Poland, 
email: atg@th.if.uj.edu.pl, jurkiewicz@th.if.uj.edu.pl 

c Institute for Theoretical Physics, Utrecht University, 
Leuvenlaan 4, NL-3584 CE Utrecht, The Netherlands, 
email: r.loll@phys.uu.nl 



Abstract 

We show that the quantum universe emerging from a nonperturbative, Lorentzian 
sum-over-geometries can be described with high accuracy by a four-dimensional 
de Sitter spacetime. By a scaling analysis involving Newton's constant, we es- 
tablish that the linear size of the quantum universes under study is in between 
17 and 28 Planck lengths. Somewhat surprisingly, the measured quantum fluc- 
tuations around the de Sitter universe in this regime are to good approximation 
still describable semiclassically. The numerical evidence presented comes from a 
regularization of quantum gravity in terms of causal dynamical triangulations. 



1 Introduction 



To show that the physical spacetime surrounding us can be derived from some 
fundamental, quantum-dynamical principle is one of the holy grails of theoretical 
physics. The fact that this goal has been eluding us for the better part of the 
last half century could be taken as an indication that we have not as yet gone 
far enough in postulating new, exotic ingredients and inventing radically new 
construction principles governing physics at the relevant, ultra-high Planckian 
energy scale. - In this letter, we add to previous evidence that such a conclusion 
may be premature. 

Our results are obtained in the framework of Lorentzian simplicial quantum 
gravity, based on the concept of causal dynamical triangulations (CDT). While 
referring to pQ EJ [3] for details, briefly, it defines a nonperturbative way of doing 
the sum over four-geometries, assembled from triangular building blocks such that 
only causal spacetime histories are included. To perform the actual summation, 
one rotates them to spacetimes of Euclidean signature. The building blocks are 
four-simplices characterized by a cut-off a, the side length of the simplices. The 
continuum limit of the regularized path integral corresponds to the limit a —>■ 0, 
possibly accompanied by a readjustment of bare coupling constants, and such that 
the physics stays invariant. The challenge of a quantum field theory of gravity is 
to find a theory which behaves in this way, and suitable observables to test it. 

How can we judge whether CDT can be taken seriously as a regularized quan- 
tum field theory of gravity? Our knowledge of the physical world suggests that a 
viable theory should generate a 'background geometry' with positive cosmological 
constant, superposed with small quantum fluctuations. The challenge is to obtain 
this from a background-independent formulation where no background spacetime 
is put in by hand. We have earlier provided indirect evidence for such a scenario 
[U [5]. Here, we present new computer simulations which confirm this picture 
much more directly, by establishing the de Sitter nature of the background space- 
time, quantifying the fluctuations around it, and setting a scale for the universes 
we are dealing with. 

2 Macroscopic de Sitter universe 

By construction, each path integral geometry (Euclideanized by the analytic con- 
tinuation mentioned above) is obtained by gluing together four-simplices in a way 
that respects a global foliation in discrete proper time t. Accordingly, each four- 
simplex will have "time-like" and space-like links of length a, our short- distance 
cutoff. Spatial slices are compact, three-dimensional manifolds of topology S 3 , 
glued from space-like tetrahedra. During simulations we fix the time direction to 
a given discrete length T to t, and for reasons of convenience take it to be periodic 
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(with sufficiently long period so as to not affect results). The bare action used is 
the so-called Regge action, a geometric realization of the Einstein-Hilbert action 
on piecewise linear geometries. 

A Monte Carlo (MC) simulation of the path integral - for computer-technical 
reasons performed at (almost) fixed four-volume - will generate a sequence of 
spacetime histories, represented by triangulations obeying causal gluing rules. 
The key result of CDT consisted in demonstrating that in a certain range of 
bare coupling constants Q this leads to the emergence of a four- dimensional uni- 
verse of well-defined temporal and spatial extension [I]. Here, we will analyze the 
geometric nature of this universe in more detail. 

Plotting the 'shape' of MC-generated path integral configurations in the form 
of a time-dependent spatial three- volume Vs(t), one typically obtains 'blob '-like 
extended geometries of time extension T ~ \/V^. Individual such configurations 
are of course not observable, in the same way that individual particle paths in 
a standard path integral are unphysical. There, one obtains the average particle 
trajectory by taking the average over paths in the path integral (with a weight 
dictated by the classical action), which agrees with the classical particle trajec- 
tory up to corrections of order h. Doing this, for instance by MC simulations, 
yields information about both the average trajectory and the size of the quantum 
fluctuations, which in general will also be of order H. 

We will follow the same procedure for our gravitational histories. Averaging 
over many different computer-generated configurations will produce the back- 
ground geometry around which the quantum universe fluctuates. Since the "cen- 
tre of gravity" of our universes performs a random walk along the time direction, 
we take the average (Vs(t)) by identifying in each history the peak of the volume 
distribution with t = 0. (T to t is always chosen much larger than the time exten- 
sion T of the universes.) The results of measuring the average discrete spatial 
size of the universe at various discrete times i are illustrated in Fig. [T] and can be 
neatly summarized by the formula 



where N 3 (i) denotes the number of three-simplices in the spatial slice at dis- 
cretized time i and A/4 the total number of four-simplices in the entire universe. 

We have verified relation (CQ) for A/4 ranging from 45.500 to 362.000. The 
overlap of curves for different A/4 according to ([1]) constitutes a beautiful example 
of finite-size scaling |6J. Eq. (TjQ) shows that spatial volumes scale according to 

1 Other choices of bare couplings can lead to different and degenerate geometries [T]. 




(1) 
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Figure 1: Background geometry (N 3 (i)): MC measurements (for fixed N 4 = 
362.000) and best fit p]) yield indistinguishable curves at given plot resolution. 
The bars indicate the average size of quantum fluctuations. 

N^ 4 and time intervals according to iV^ 4 , as one would expect for a genuinely 
/bur-dimensional spacetime. Translating ([I]) to a continuum notation leads to 

V^3 d (*) = ^cos 3 (|), (2) 

where we have made the identifications 

V 4 = a 4 C 4 N 4 , Jg7 t V 3 = a 3 C 4 N 3 , U = a i. (3) 

In d2J), C4 = \/5/96 from the discrete four-simplex volume [2], and y/g~^ is the 
constant proportionality factor between the time t and genuine continuum proper 
time r, t = ^fgTt t. Writing V4 = 87r 2 i? 4 /3, and ^fg~tt = R/B, eq. ([2]) is seen to 
describe a Euclidean de Sitter universe (a four-sphere, the maximally symmetric 
space for positive cosmological constant) as our searched-for, dynamically gen- 
erated background geometry! In the parametrization of ([2]) this is the classical 
solution to the action 

+ (4) 
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where k 2 = 9(27r 2 ) 2 / 3 and I is a Lagrange multiplier, fixed by requiring that 
the total four-volume be V4, / dt^Jg~^ t V 3 {t) = V4. Up to an overall sign, this is 
precisely the Einstein-Hilbert action for the scale factor a(t) of a homogeneous, 
isotropic universe (rewritten in terms of V 3 (t) = 27r 2 a(t) 3 ), although we of course 
never put any such simplifying symmetry assumptions into the CDT model. The 
discretized, dimensionless version of (HI) is 



(N 3 (i + 1) - N 3 (i)) 2 ^ ~ h ATl/3 
N 3 (t 

where k 2 oc k 2 . The identifications ([3]) lead to a naive continuum limit of the 
discretized action S 'discr with 

a 2 C 

Our next aim will be to determine the coefficient k\ in front of the effective 
action (JSJ) from the computer simulations. Because of relation ([6]), this will give 
us an estimate of the gravitational coupling constant G in terms of the lattice 
spacing a. Since in our units the Planck length is £pi = \/G, this will set a 
physical length scale for lattice quantities. While the classical solution (j2J) does 
not provide information about the numerical value of k\ in front of ([11) , a saddle 
point calculation shows that the fluctuations SV 3 (t) := V 3 (t) — V 3 (t) around 
V 3 d (t) = (V 3 (t)) will be of the order 

((8V 3 ) 2 )~GV 4 . (7) 

Therefore, if it is true that also the fluctuations SN 3 (i) := N 3 (i) — N 3 (i) of 
our model are well described by the mini-superspace action (jlj), we can simply 
determine k\ from measuring their correlator. 



3 Fluctuations around de Sitter space 

Having demonstrated that the action (j3J) gives a perfect description of the mea- 
sured V 3 cl (t), we will now show that it also describes the observed quantum fluc- 
tuations around de Sitter space, defined by the correlator 

C{t,t') = (8V 3 {t)8V 3 {t')). (8) 

Our first observation from the data is that the discretized version of ([H]) scales 
with the four-volume according to 

C N4 (i,i') = (5N 3 (i)6N 3 (i')) = N 4 F^/Nl'^i/Nl 1 ^ (9) 
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Figure 2: Analyzing the quantum fluctuations of Fig.[TJ diagonal entries F(t, t) 1 / 2 
of the scaling function F from ©, for iV 4 = 45.500, 91.000, 181.000 and 362.000. 

where F is a universal scaling function of order a . This is illustrated by Fig. [2] for 
C^ 2 (z, i), corresponding precisely to the fluctuations ({8V 3 {t)) 2 ) 1 / 2 of Fig. [IJ This 
scaling implies that the action ([5]) can only describe the fluctuations measured 
for different N± if k\ is independent of N±, thus confirming the scaling behaviour 
G ~ a 2 anticipated in ([6]). 

To demonstrate that F(t,t') is indeed described by the effective actions (J4]), 
([5]), let us for convenience adopt a continuum language and compute its expected 
behaviour. Expanding (jlj) around the classical solution as V 3 (t) = V£ l (t) +x(t), 
the quadratic fluctuations are given by 

W(Mt <)> = M e-U/™...>,'> 

= M-\t,H), (10) 

where T>x(s) is the normalized measure and the quadratic form M(t, t') is deter- 
mined by expanding the effective action S to second order in x(t), 

S(V 3 ) = S(V 3 d ) + I dt x(t)Hx(t). (11) 
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In (11 ip . -H" denotes the Hermitian operator 



H dtcos 3 (t/5)dt 5 2 cos 5 (t/5)' ^ 

which must be diagonahzed under the constraint that J dtJlfrt x(t) = 0, since V4 
is kept constant. 

Let e n (t) be the eigenfunctions of the quadratic form given by ( TTTT) with the 
volume constraint enforced, ordered according to increasing eigenvalues l n . If 
this cosmological continuum model were to give the correct description of the 
computer-generated universe, the matrix 

M-\t,f) = J2 6n{t ] en{t ' ) . (13) 

n=l 

should be proportional to the measured correlator C(t,t'). Fig. [3] shows the 
highest eigenfunction calculated from the data, the matrix C(t,t'), and the cor- 
responding lowest eigenfunction calculated from the effective action, the matrix 
M(t,t'). The agreement is very good, in particular taking into account that no 
parameter is adjusted in the action (we simply take B = sqN^ 4 in (|2J) and (jTTI) . 
i.e. 14.47a for N4 = 362.000). One can also compare the data and the matrix 
M _1 (£,£') calculated from (fT3|) directly. The agreement is again good, although 
less spectacular than in Fig. [31 Awaiting publication of a full analysis of the 
data [6], suffice it to say that within measuring accuracy the fluctuations of the 
spatial volume are described by the mini-superspace action (111) . This enables us 
to numerically estimate k\ ~ 0.016, which according to ([6]) leads to G ~ 0.22a 2 , 
or £pi « 0.47a. In other words, the linear size of the quantum de Sitter universes 
studied here lies in the range of 17-28 Planck lengths. 



4 Discussion 

The CDT model of quantum gravity is extremely simple, namely, the path integral 
over the class of causal geometries with a global time foliation. In order to perform 
this summation explicitly, we introduce a grid of piecewise linear geometries, much 
in the same way as when defining the path integral in quantum mechanics. Next, 
we rotate each of these geometries to Euclidean signature and use as bare action 
the Einstein-Hilbert action H in Regge form. Nothing else is put in. 

The resulting superposition exhibits scaling behaviour as function of the four- 
volume, and we observe the appearance of a well-defined average geometry, that 

2 Of course, the full, effective action, including measure contributions, will contain all higher- 
derivative terms. 
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Figure 3: Comparing the highest eigenvector of C(t, t') and the lowest eigenvector 
of M-'(U'). 

of de Sitter space. We are definitely in a quantum regime, since the fluctuations 
around de Sitter space are sizeable, as can be seen in Fig. [TJ Both the average 
geometry and the quantum fluctuations are well described by the mini-superspace 
action Unlike in standard cosmological treatments, this description is the out- 
come of a nonperturbative evaluation of the full path integral, with everything 
but the scale factor (equivalently, Vs(t)) summed over. Measuring the correla- 
tions of the fluctuations in the computer simulation enabled us to determine the 
continuum gravitational coupling constant to G ~ 0.22a 2 , thereby introducing 
an absolute physical length scale. Within measuring accuracy, our de Sitter uni- 
verses (with volumes in the range 22.000-173.000 £pi) are seen to behave perfectly 
semiclassically. 

Can we study smaller universes, which are themselves of Planck size? Taking 
the coupling G as a true measure of the gravitational coupling constant, the 
simplest way is as follows. We are free to vary A/4 and the bare gravitational 
coupling constant g of the Regge action (see [T] for further details on the bare 
coupling constants), with the effective constant k± a function of go- If we adjusted 
go such that in the limit — > oo both 

V A ~N A a A and G ~ a 2 /kx(g ) (14) 

remained constant (i.e. k\(go) ~ l/x/AQ, we would eventually penetrate into the 



8 



(sub-)Planckian regime. Since we have already seen deviations from classicality 
at short scales [JJ, we would expect the canonical scaling of G to change there, 
or, stated differently, the simple effective action (j3J) to be no longer valid. Renor- 
malization group methods have produced predictions for the scaling violations of 
G in the context of asymptotic safety [7] , which in principle we should be able to 
test. In this context it would be ideal to have an observable with an associated 
correlation length that could be kept constant when expressed in terms of V^ A . 
A further step will be to include matter in the model and verify directly that G 
can indeed be interpreted as Newton's constant, perhaps along the lines pursued 
earlier in Euclidean quantum gravity [8]. All of these issues are currently under 
investigation. 
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